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I. INTRODUCTION 



■ Renormalization Group (RG) flow equations pioneered by Wilson and Kadanoff allow for a prediction on the 
behavior of a interacting theory in different momentum regimes [Q. When probing physics in the infrared (IR) it is 
' often useful to derive a low-energy effective theory by integrating over the irrelevant short-distance (high momenta) 
^ ' modes. Thus one is left with a so-called Wilsonian or low-energy effective action distinguishing it from the one- 

particle-irreducible (IPI) effective action. 
\l In order to separate the fast-fluctuating short-distance modes from the slowly- varying ones one has to introduce 

^ a momentum scale fc, which appears naturally in the momentum cutoff rcgularization. The momentum cutoff reg- 
ularization can be formulated by means of a blocking transformation. This procedure is achieved by introducing a 
' so-called smearing or IR cutoff function fk, which suppresses the fast-fluctuating modes of the original unblocked 
fields In that way the fields are naturally separated into fast and slow mode components. The integration over the 
fast modes corresponds to the blocking transformation on the lattice. This results in an effective action parameterized 
by the averaged blocked field at the scale k. The full IPI Feynman graphs or vertex functions, providing the natural 



I ■ tool for the study of broken symmetries in interacting field theories, are generated by the renormalized effective action 
Oh in the limit A: — *■ 0. Thus the /c-dependent effective action provides a smooth interpolation between the bare action 
^ ' defined at the UV scale A, where no quantum fiuctuations are considered and the renormalized effective action in 
the IR. The RG fiow pattern of the theory is obtained by studying an infinitesimal change of the IR scale k, in the 
effective fc-dependent action. A nice introduction to this topic can be found in the recent lecture 

Unfortunately, this integration step cannot be computed in an exact way and one has to perform some approxi- 
mation like the frequently used loop expansion. Based on an one-loop expression with a heat-kernel (proper-time) 
' rcgularization we derive renormalization group improved flow equation and study different expansion pattern. 

It is well known that a sharp momentum cutoff rcgularization is in conflict with the gauge invariance [^j^. The 
crucial task is to implement both the UV and the IR cutoff scales without destroying the important symmetries. For 
such reasons we use the operator cutoff rcgularization in Schwinger's proper-time representation ||^ which neither 
depends on the dimension nor violates any symmetry of the theory. 

The paper is organized as follows: in Section ^ we define the model Lagrangian and give a general introduction 
to the renormalization fiow equation using the heat-kernel rcgularization and different choices of the cutoff function. 



In Section [II we approximate the potential by a polynomial expansion and derive the zero temperature coupled 



fiow equations. Here we also present the diagrammatic interpretation of the equations. Temperature is introduced 



in Section IV and the fiow equations associated with the expansion of the potential are shown for different cutoff 
functions. In Section ^ we present our results at zero temperature, compa ring the results of the expansion technique 



at different orders to exact grid solutions of the fiow equations. In Section the finite temperature calculations are 



presented and we discuss the sensitivity of the critical temperature on the expansion parameter, while in Section VII 
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the behavior of the critical exponent, /3, is studied. Finally, in Section VIIl we conclude by critically assessing the 
feasibility of finite-density calculations using the expansion technique. 



II. RENORMALIZATION GROUP FLOW EQUATIONS 

Since the intention is to study the dynamical breaking of chiral symmetry by matching the short-distance regime of 
QCD with the large-distance behavior in the framework of evolution equations we consider an effective theory which 
incorporates simultaneously quark fields and composite mesonic fields as the active degrees of freedom. A suitable 
choice is the linear sigma model with a- and 7r-mesons coupled to Nf = 2 constituent quarks. The partition function 
in the chiral limit and without external sources is given by the Euclidean-space path integral 

Z = J-DqVqVaVnexp{- Jd^x {£f + Cb)}, (1) 

with a fermionic part including a Yukawa coupling constant g, 

Cf ^qhd + g{(J + ifTrj5))q, (2) 



and a bosonic Lagrangian part 



Cb = 1 {id,a)^ + {d,nf) + V{a' + 7?^) , (3) 



respectively. The mesonic self-interaction is parameterized by the potential term V. At finite temperature the integra- 
tion in the action along the temporal direction is understood on the stripe [0, P] where /3 is the inverse temperature. 
The parameters of the linear cr-model are fixed at zero temperature as in ref. |@]. Because the theory is strongly 
coupled we cannot calculate the low-energy theory directly. We assume that at an ultraviolet scale Awl GeV, the 
full QCD dynamics reduces to a hybrid description in terms of (massless) current quarks and chiral bound states. 
The gluons are frozen in the residual mesonic degrees of freedom, their couplings and the effective potential V\ which, 
at this scale A, takes the following form: 

V^{^) = '^$' + ^$\ where $={cT,n) (4) 

is an 0(4)-symmetric vector. Both and the four-boson coupling A, are positive. During the course of iteratively 
including quantum fluctuation from a restricted scale range the potential and accordingly all couplings become k- 
dependent. In this work we neglect the evolution of the Yukawa coupling g, and set it equal a constant value of 3.2 
in order to reproduce the constituent quark mass Mq « 300 MeV, with the pion decay constant = 91 MeV, in the 
infrared. 

The one-loop contribution to the effective action yields in general a non-local logarithm, which is ultraviolet divergent 
and needs to be regularized. This regularization is based on the general proper-time representation introduced by 
Schwinger . After introducing a complete set of plane waves we are lead to the following one-loop correction to the 
effective action for the bosons 



1 f, 7dT r d^q --U'+^-^9+d'+j0^^ 

d X — / , , Tr e 



2 J J T J (27r) 



(5) 



and similarly for the fermions 







where we have split the linked integration over fermions within the bosonic integration into two parts, evaluated at 
the minimum of the effective potential. Thus the total effective action becomes a sum of both terms F = F^ + F^ 
each depending on 0o- For details we refer the reader to ref. 0. 

A prime on the potential denotes differentiation with respect to the fields 0^ and we elucidate our notation by the 
equation 
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Tr e"" W7 = Se^^-^' + ^-r{2V'+A4.-v") 

The short-distance or ultraviolet divergences appear as divergences at r = and the large-distance or infrared 
divergences as divergences at r ^ cxd. 

In spite of a continuous blocking transformation procedure we modify the above expressions by introducing a 
regulating function fk (r) into the proper-time integrand which separates the "fast" and "slow" fluctuating modes as 
well as respects the important symmetries of the considered theory Of course, the form of the regularized flow 
equations does depend on the choice of the a priori unknown blocking function fk- In order to find an explicit form 
for fk one is lead by the following conditions: 

(1) the modified fc-dependent effective action Tk should tend to the full generating functional F in the IR limit 
fc — > 0, i. e. when the infrared cutoff /fc, is removed. Thus we require that fk^o^r —>■ 00) — > 1. As a consequence in 
the limit fc ^ all quantum fluctuations are taken into account. 

(2) in addition to the previous condition which regularizes the IR we set fk{T = 0) = l for arbitrary k. This means 
that the UV regime in the proper-time (r = 0) is not regularized by this requirement. Since we start the evolution at 
a finite (large) UV scale A, an UV regularization through the cutoff function is not necessary. 

(3) a differential equation for fk can be deduced by comparing the heat-kernel regularization to a momentum cutoff 
regulator in the one-loop contribution where the one-loop trace is restricted over a momentum slice k < p < A. It 
readsQ: 

fk ^ -{rk^ fhirk') (8) 

with h{x) being any regular function around the origin (cf. 0,^). Consider for example the RG improved leading-order 
contribution to the one-loop potential of a scalar one-component theory in a sharp cutoff regularization, within the 
momentum shell k < p < A, (for simplicity we show it for the second term of Eq. (^, the first term can be treated 
similarly.) 

A 

^'^IJ (§4ln(p' + 2U' + 40V") . (9) 

k 

Differentiating this equation with respect to k results in the flow equation 



dk (47r) 



k^ r dr 



(47r)2 J T 







where in the last line we have inserted the proper-time representation of the logarithm.^ 

This expression can be compared to the modified heat-kernel regularization and yields the required differential 
equation of the type (||) 

kf'k = -2(Tfc2)2e-fe^ . (11) 
It is also possible to generalize this result to arbitrary d dimensions with the result 

kfk - -^i^kY/^Y^jYf'^'^ ■ (^^^ 

One suitable solution of Eq. ( pi] ) which fulfills the proper boundary conditions is the function 

fi'\Tk^)^e-^'^\l+rk^) . (13) 



prime on fk means differentiation with respect to k. 
^Proper regularization is assumed but not indicated in this formula, see e.g. [|5l|l 
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We note, that despite the evohition equation derived from a sharp momentum cutoff and the one from the heat-kernel 
regularization with the smooth cutoff function (|l3| ), have the same form, this neither means that the two regularization 
methods are equivalent nor that the choice (|13|) is unique. E.g. this choice for fk is different from that in ref. 0. 
In fact, it provides a slower IR convergence as fc — » 0. Including higher monomials of the form rfc^ in the cutoff 
function fk will accelerate the IR convergence The choice of the cutoff function affects the explicit form of the 
flow equation itself. To show this we present our result with both the cutoff ( |l3|) (I), and the one with higher orders, 

n 

/(")(rfc^) = e-'=^^i(rfc^r. (14) 

4=0 

where n = 2 corresponds to the original cutoff of ref. [Q fj,^^^ . We shall also examine the cutoff f^^^^^ corresponding 
to n = 3. The introduction of these higher-order cutoff functions turns out to be even necessary in the case of higher- 
order gradient expansion in order to attain equality between the momentum cutoff and the proper-time regularization 
on the level of the wave function renormalization [^^. 

Differentiating the effective action or potential Vk which become /c-dependent through the introduction of the 

(n) 

blocking function ff, with respect to the arbitrary infrared scale k and after calculating the trace over inner spaces 
in Eqs. (^J^) yields flow equations which incorporate only fluctuations from one- loop order and will break down in 
the high temperature regime pT[ . 

This so-called "independent-mode approximation" can be further improved by taking higher graphs like daisy and 
superdaisy diagrams into account. This can be accomplished by considering the interactions between the fast 
and slow modes and consists of replacing the bare potential V on the rhs of the flow equation with the fc-dependent 
potential Vk (cf. e. g. Eq. (^)). 

Taking this substitution of the full potential Vk into account we find the following flow equations for the potential 
in 4-dimcnsion (cf. with ref. [l^ ) 

- -j^ {31n(l + 2y,7fc2) + Hl + 2V^/k^ +M'Vl'/k^) 

-ANfNAn{l+9^4>Vk^)} , (15) 

for cutoff function (I) and 

dv!:"'> 



k- 



dk n(n-l)167r2 [(l-H2F;!/A:2)«-i (1 + 2^/^7^ +402v;"//c2)r 
47V/iVc 

~(l-|-g202/fc2)n- 



(16) 



for the other cutoff functions (n = 2 corresponds to case (II) and n = 3 to case (HI)), respectively. Using the smooth 
cutoff function for case (I) /^.^^ yields a Renormalization Group equation with a characteristic logarithmic structure 



very similar to the Wegner-Houghton equation with a sharp-cutoff |14|. 

The first term on the rhs of Eqs. (p"5|-[l6|) represents the flow equation contribution from the three massless pions^, 
the middle term the sigma-meson contribution with mass squared term 2y^'+4(/)2y^" and the last term the NfNc quark 
contribution]^ 

As a result of the choice of the cutoff function /^."'' for n > 1 (cf. ref. |0]) the flow equation contributions have 
the typical form l/{k^ + mass2) and are sometimes called "threshold functions" |l^. This structure controls the 
decoupling of the massive modes in the infrared. In the case of cutoff function (I) however, the "threshold functions" 
are represented by logarithms. Nevertheless, the appearance of the mass terms serve here the same purpose, namely 
to regularize the evolution equations. 

At first glance the structure of the nonperturbative flow Eqs. (p^p^ looks different. Eq. (^ is the one-loop 
resummation of the perturbative expansion and can be rewritten by expanding the logarithm in the non-Gaussian 
pieces of the potential. Comparing the resulting series with the expansion of the other flow Eqs. (^6|) one finds 
similar terms relating Eqs. (Eq-Qq) to each other. 



^14' vanishes at the minimum 

^particle - antiparticle and spin 1/2 yields the factor 4. 
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On the other hand one can expand the effective action in powers of momenta. In coordinate space this expansion 
takes the form (gradient expansion) 

r[0o] = J d^x |-y[0o] + iz[0o](9^</'o)' + r[0o](9,,0o)" + • • j ■ (17) 

The lowest-order approximation of this expansion is the local potential approximation (LPA) which consists of con- 
sidering a constant background field, setting Z = 1 and neglecting all higher coefficients functions Y^ --. Now the 
truncation becomes apparent and one sees the restriction to the effective potential neglecting the influence of the 
anomalous dimension -q defined by —k ^^Qif - A more general investigation of an 0(iV)-symmetric potential including 
anomalous dimensions within this approach will be published elsewhere pi| . 



III. EXPANSION OF THE POTENTIAL 



In the following we introduce a short-hand notation p = 0^, and expand the potential Vk at the scale k around its 
minimum to study the stability of such an expansion as the function of order M . A similar analysis involving only 
mesonic fields was done recently in ||2^ . 

Explicitly, for the symmetric phase (po = 0) up to order M, 



^ 1 



^t'^-Y.7^^-^r.{k)p- (18) 
and similarly for the broken phase around the new fc-dependent nontrivial vacuum (po 7^ 0), 

M 

= ^" + E (^^2„(fc) (P-Po(fc))" . (19) 

The second coefficient &2 of this expansion defines the minimum po of the potential and therefore vanishes. In the 
local potential approximation the expansion coefficients of the potential correspond to the 2n-point proper vertices 
evaluated at zero momenta which we have denoted by a2n for the symmetric phase. Now substituting the upper 
expansions on both sides of the flow equation disl ) we can deduce a coupled system of flow equations for the proper 
vertices. For the symmetric phase (evaluated at po = 0) we obtain the following set of coupled equations for the first 
four couplings ag, a2, 04 and ag with the cutoff function (I), 

^'-^ = ^.'^T^ (20) 
dk 47r2 fc2+a2 

k^ f 10 as a^oe 100 (4__ QanMArS^X ^o^-a 

Note, that the flow equation of the first coupling oq, has a contribution from the quarks in the symmetric phase 
through a2 ■ Similar equations may be obtained using the cutoff functions (II) and (III) . 

The integration of the high-momenta degrees of freedom results in new one-loop flow contributions to the couplings 
of the effective Wilsonian Lagrangian. During the evolution towards the IR all possible higher-dimension composite 
interactions are generated which are even partially nonrenormalizable. However, the contributions of these non- 
renormalizable vertices are under control, since the scales involved during the evolution are much smaller than the 
corresponding cutoff values. Hence, these operators have a negligible effect on the physics at infrared scales. 

Equations (pl|-|2^) are easily identified as the differential form of the truncated Dyson-Schwinger equations at 
vanishing external momenta. The denominator + 02 plays the role of the bosonic propagator with the boson self- 



energy squared 02 (fc) and Eq. (21) corresponds to the bosonic self-energy equation at vanishing external momenta. 
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(24) 



^« -M 




while Eq. ( p^ ) is the Bethe-Salpeter equation for the vertex function at vanishing external momenta, 

-:■< - ;<:::< - ( ) 

without the integration in the loop-momentum variable k. In the figure the single external lines represents the 
low-momenta field (which do not appear in the equations), while the inner lines are dressed by the high- momenta 
fluctuations in one-loop corrections down to scale k. Each of Eqs. (pT|-p3|) symbolizes the sum of all one-loop diagrams 
with a fixed number 2n of external legs |po| . The analogous fermionic inverse propagator is associated with a massless 
quark in the symmetric phase, and is proportional to k. 

For the broken phase we obtain the following flow equation structure expanding around the nontrivial minimum 
Po, 

= + ^4Po/3fc2) - iN^Nf log(l + g'po/k')} (26) 



dk (47r)2 

(Gbi 2 
+ 5 

dbi _ k^ (3b6 hi 1 (1564 + Pobe? 1 3566 + Pobs 



jdpQ_ fc4 f 664 2 1564 + po66 48g2jVe7V/ ; 

dk (47r)264 [k^ 5k^ + 64P0/3 P + g2p^ r y ) 

k 



dk (47r)2 {5k^ k^ 75 (fc2 + pab^/if 35 (fc2 -|- P064/3) 



[k^ + g^Pof] ' 10 dk ' ^ ' 

We omit here the explicit form of the higher flow equations because the result is lengthy but straightforward. If one 
defines the mass of the cr-particle as m'^ = po&4/3 then one can easily identify the massive cr-propagators which must 
be distinguished from the massless pion propagators in the chiral limit. On the other hand the quarks are now massive 
m2 = g2pg in the broken phase. Furthermore one sees nicely the matching of both equations (|2|) and ( |2^ ) during the 
course of evolution with respect to the scale k because, at the chiral symmetry breaking scale, 02 and po tend both 
to zero. This feature also applies to the flow equations for higher couplings be etc.. 
Once again the diagrammatic interpretation of Eq. ( ^ ) is straightforward, 



dk 



oc 



- <-> + C ) 



(29) 



The crosses on the graphs symbolize the vacuum expectation value po- The massive cr-propagator is represented by a 
double line, the massive 7r-propagator as a dashed line and the massive quark-propagator as a thick solid line. Below 
several graphs the order of the corresponding couplings is indicated. 

We note, that the form of the higher-order cutoff functions in the heat-kernel regularization corresponds to 
adding n — 1 ghost fields H to the loop integrals i.e. additional graphs appear in Eqs. (|24|j25|,|29|) . 
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IV. FINITE-TEMPERATURE FLOW EQUATION 



It is straightforward to generalize the above approach to finite temperature. In order to match the four-dimensional 
theory at zero temperature with the effectively purely three-dimensional behavior at the critical temperature it 
is more convenient to use the imaginary time approach (Matsubara formalism). This means that we replace the 
four-dimensional momentum integration in Eqs. (5|-^) by a three-dimensional integration and introduce a discrete 
Matsubara summation with the corresponding frequencies for bosons w„ = 2n7rT, and fermions Vn = (2n+l)7rr, in 
the zero- momentum component. 

Furthermore, an additional length scale /3, the inverse temperature, is introduced. This incorporates additional 
thermal fluctuations beyond the already existing quantum ones. As shown in the appendix of ref. it is possible to 
split the thermal fluctuations into a zero temperature contribution plus a finite temperature contribution using the 
so-called generalized 0-function transformation. However, in the present work we decided to solve the Matsubara 
sums numerically and do not apply the 0-function transformation. 

Thus at flnite temperature the scale k serves as a generalized IR cutoff for a combination of the three-dimensional 
momenta and Matsubara frequencies. A careful investigation of how the two types of modes are controlled by different 
choices of the blocking function can be found in ref. p^ . In the present work we use a cutoff function family with a 
four-dimensional momentum variable, related to the smeared version of type (3) in ref. p3t . 

Using the cutoff function (I) - Eq. ([l^) - we obtain the following flnite-temperature (dimensional) flow equation 
for the full potential 

ay(^) p ~ r 3 1 

^ 8^^„£o I [l + (^^ + 2V;0/fc2]i/2 + [l + K2+2V;!+402F,")/fc2]i/2 

(30) 



AN.Nf 



and 



[l + (z/2+5202)/fc2]l/2 



jdvt^ k\ (2»-3)!! 

dk Stt n! 2"~i ^f-'^ \ [l + (u;2-H2Vfc')/fc2]"-i/2 



1 m^Nf 



[l + (L^2+2V^fe' + 402V;")/fc']"-'/' [l + (i/2+g2^2)/fc2]«-l/2 



(31) 



with (2n— 3)!! — 1-3 . . .•(2r7r- 3) for cutoff II {n — 2) and III (n — 3). Due to the three-dimensional momentum integration 
fractional powers arise in the threshold functions. In appendix A of ref. a relation for the low-temperature limit 
of the Matsubara sums to the four-dimensional integrals can be found. This guarantees the correct matching of the 
finite-temperature equation to the zero temperature one since in the limit T — > the two become identical. 



V. RESULTS AT T = 0. 

At zero temperature the general flow equations (30-^l|) reduce to the flow equations (|l5|-[l^) which can be proved 
analytically using the low-temperature limit of the Matsubara sums as described in the appendix of ref. Q . 

We have solved numerically the flow equations (^0|-|3l|) with the cutoff functions (I)- (III) both on a grid and using 
the set of coupled equations on the expansion coefficients (p0|-p3|) and ( p6|-p^ ) and compared the results. In the grid 
calculation the potential was discretized by 80 points in the range < (jrK 0.05 GeV^ and the set of 80 coupled 
differential equations were solved with a Runge Kutta method. 

We start the evolution in the symmetric phase deep in the ultraviolet region of QCD at A 1.2 GcV and fix our 
parameters similarly to ref. 0| for cutoff function (II) , to get at the end of the evolution the physical value of the pion 
decay constant, and the chiral symmetry breaking scale in the range of 0.8 — 1 GeV. 

In table | we list the used parameterization for the cutoff functions (I), (II) and (III). In the course of the evolution 
towards zero with respect to k we encounter the transition to chiral symmetry breaking at a finite scale k^sB, where 
we switch to the equations for the broken phase. The chiral symmetry breaking scale depends on the underlying 
cutoff function and is also listed in table |. The required initial values are chosen in such a way as to fix the pion 
decay constant /^r, at the end of the evolution in the infrared. 

In Fig. 1^ we compare the convergence of the expanded equations as a function of the order of the expansion. 
The solid lines represent the grid value of /jr and the triangles are the ones related to the solution of the expanded 
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equations for the cutoff I (left), II (middle) and III (right). A similar plot for the chiral symmetry breaking scale 
kxSB is presented in Fig. ||. 

The lesson we learn from the figure is that for the cutoff (I) the convergence is very poor, actually the higher-order 
calculations are getting even worse, while for the higher-order cutoffs the convergence is improving rapidly with the 
order of the cutoff n. The lowest-order M — 2, truncation overestimates /^r by 70% for cutoff (I), by 25% for cutoff 
(II) and only by 8% for cutoff (III). Cutoff (I) shows unstable behavior when increasing the order of expansion, while 
cutoff (II) gives already a good approximation to the full (not expanded) calculation at M = 6 and cutoff (III) at 
M = 4. However, for all cutoffs, taking higher-order polynomials, the global shape of the potential becomes worse 
going away from the minimum, that is the convergence radius of the approximation becomes smaller and smaller 
around the minimum. This is demonstrated by calculating the infrared limit of the expansion coefficients. 

As mentioned above, we fix the parameters of the quartic potential (^) at the momentum scale A — 1.2 GeV, setting 
the initial conditions for the evolution equations as oq = 0, 02 = m^, 04 = 6A and a2i — for z > 3. At the chiral 
symmetry breaking scale k^sB, the mass parameter 02, becomes negative and the minimum of the potential shifts to 
some finite po value. Below this momentum scale we use Eqs. (p6|-p8|) for the broken phase. Studying the structure 
of these equations and higher orders in the k ^ limit, one finds 



'dk 



36 bl 



yielding & 



{III) 



36 In A: 



and 



Generally, 



566 
dk 



768 hi 



yielding 6, 



(///) 



2^4 



1215fc2 1n^yfc 



1 



(32) 



(33) 



(34) 



This result is also supported by the numerics. Since the asymptotic structure is similar for all cutoff functions, we get 
similar results. 



1^ and 6^-)--!^ 



Infc 



In k 



(35) 



for the cutoff I and II. 

The expansion coefficients for 2i > 6 diverge in the infrared limit, hence making the expansion of the potential 
meaningless outside the minimum. 



VI. RESULTS AT T / 0. 

In Fig. 1^ we present the critical temperature as the function of the order of polynomial expansion for the same 
initial parameter set as for zero temperature. Once again there is no convergence obtained for cutoff (I), however, 
cutoffs (II) and (III) behave properly. The grid value of the critical temperature is Tc « 150 MeV for all three cutoffs. 
For cutoff (II), M = 2 overshoots the asymptotic value by 25% and stabilizes at M > 6. Cutoff (III) is much more 
stable, deviating only by 9% from the grid value at Af = 2 and stabilizing at M = 4. 

The evolution of the minimum of the potential at the critical temperature shows the same interesting and coun- 
terintuitive pattern already discussed in 0: decreasing the momentum scale starting from the symmetric phase in 
the deep ultraviolet one enters the broken phase due to the presence of the quarks and the minimum of the potential 
4'k = firik), first increases up to fc « 0.5 GcV, then decreases back to zero (see left side of Fig. Such a behavior 
is confirmed by the calculations made on the grid, the shape of as a function of k can well be fitted by an ellipsis. 
Above the critical temperature lowering the momentum scale hence there are two transitions, one from the initial 
symmetric phase to a broken one at k^sB^ and at small k another transition is driving the system back to the sym- 
metric phase. This indicates that the infrared modes now restore the symmetry. Such a behavior was also observed 
in other works with a different cutoff family [T^j2^ . 

There is an interesting relation that can be found between the value of the pion decay constant at zero temperature 
and the critical temperature for the same initial parameter set but different expansion orders. The obtained values 
follow nicely a linear fit (see left part of Fig. ^), in accordance with a slope parameter of 2 found in |^ in the case 
of the NJL-model. 
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The asymptotic behavior of the expansion coefficients b2i, may be studied analytically in the fc ^ limit of the 
broken phase. Since in this limit the relevant quantity T/k is approaching infinity, the Matsubara sum may be replaced 
by its lowest-order contribution. Keeping the leading-order terms in the momentum scale k, we arrive at 

for i > 1 and up to T < Tc. Once again the coefficients are divergent for i > 3 shrinking the convergence radius of 
the expansion to one point. 



VII. CRITICAL EXPONENTS 



Our numerical results show a second-order phase transition in the order parameter , that is the order parameter 
vanishes continuously with a critical exponent /3, /^r {Tc — T)'^ . A typical plot, for M = 8, is presented in Fig. ^ 
(right), in the scaling window of the model. The critical exponent /3, is only slighlty sensitive to the order of the 
polynomial expansion and cutoff used and agrees with the three-dimensional scalar 0(4)-Heisenberg model value, 
f3 ~ 0.4. For comparison lattice Monte Carlo simulations yield B = 0.3836(46) 4 — e expansion (3 — 0.38(1) 
[0 and the average action approach f3 = 0.407 In fig. || we demonstrate the order dependence M of the 

critical exponent for the different cutoff functions (I)- (III) and compared our values with the one in the literature: the 
shadowed region represents the value region obtained in different calculations. Despite the fact that the convergence 
of the pion decay constant and critical temperature was poor (even absent) in the case of the cutoff function (I), the 
critical exponent /3, is still quite well in agreement with the 0(4)-model value. For the higher cutoff functions (II) 
and (III) the convergence becomes more and more stable, however, a systematic decrease with the order n, of the 
cutoff functions is observed in the direction of the value obtained by the 4 — e expansion and the MC calculations. 
This decrease persists further at M = 4, n = 4 approaching closer to the 4 — e expansion value. 

It is also interesting to compare the size of the scaling window for different cutoff functions. At Af = 4 we found that 
the window starts at the values of (T^ - T)/Tc = 17%, 7%, 2.7% and 0.8% for cutoffs (I) through (IV), respectively. 
This indicates a widening of the scaling regime for lower monomial cutoff functions. 



VIII. CONCLUSIONS 



We have studied a model Lagrangian of quarks and mesons with renormalization group equations based on an 
one-loop expansion within the local potential approximation, for different cutoff functions. Expanding the meson 
potential in polynomials of degree M around its minimum, we found that the lowest order cutoff function (I), 

related to the Wegner-Houghton equation with a sharp cutoff has serious convergence problems. Increasing the order 
of the cutoff function applied, the convergence in the order of the potential expansion M, increases rapidly, already 
with the next cutoff function (II), the expansion in the potential is within 25% for M = 2 {(f)'^ theory). In order to 
get results closer to the full potential (grid) calculation however, one has to take higher corrections into account. The 
analysis of the critical behavior in ref. has shown that the "magnetization" (j)o, is related to the external source 
J, (the current quark mass in the present model) as (f>o ~ J^^^ with the critical index S « 4.8. Hence the effective 
potential has the dependence Vefj ^ 0^"'"'' around the minimum. This explains the necessity to have at least terms 
with M = 3 corresponding to 0^ in the effective potential near Tc- Our study shows that for cutoff function (II) with 
M = 6 expansion yields stable results and close to the full grid simulation for several global quantities (such as the 
pion decay constant, the chiral symmetry breaking scale k^sB, and the critical temperature Tc). However, for local 
quantities, such as the potential itself, increasing order gives less and less reliable shape (the convergence radius of 
the expansion defines a smaller and smaller interval around the minimum), deviating from the full grid calculation. 
This means, that studies involving first-order phase transition (as encountered at finite baryon density [p6|j2^ ]) cannot 
be treated in such a model up to the infrared scale, since the two minima are separated further than the convergence 
radius and cannot be calculated reliably. The problem may be addressed using full grid calculations or other expansion 
basis may be considered, however, for finite density such a basis should include non-analytical functions. 
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FIG. 1. Dependence of the pion decay constant /,r, on the order of expansion for different cutoff functions. The solid lines 
represent the grid values. 
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FIG. 2. Dependence of the chiral symmetry breaking scale k^SB, on the order of expansion. The solid lines represent the 
grid values. 
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FIG. 3. Dependence of the critical temperature on the order of expansion for the three cutoffs. The solid lines represent the 
grid values. Note the cut on the y-cixis. See text. 
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TABLE I. The initial parameters rn and A at the scale k = A for different cutoff functions (1) (111) and the resulting pion 
decay cou-staut /„. at the sc-alo A- — 0. The cliiral .syimuetry hreaJduf? scale A\.s-c. is also giveu. 



cutoff function 


m [GeV] 


A 


A [MeV] 


fexss [GeV] 


(I) 


0.4 


.51.5 


91.85 


0.87 


w 


0.1 


;i().o 


l)i.8() 


0.98 


(III) 


0.4 


16.0 


91.86 


1.00 
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